Practality feature in oil price fluctuations 
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The scaling properties of oil price fluctuations are described as a non-stationary stochastic process 
realized by a time series of finite length. An original model is used to extract the scaling exponent 
of the fluctuation functions within a non- stationary process formulation. It is shown that, when 
returns are measured over intervals less than 10 days, the Probability Density Functions (PDFs) 
exhibit self-similarity and monoscaling, in contrast to the multifractal behavior of the PDFs at 
macro-scales (typically larger than one month) . We find that the time evolution of the distributions 
are well fitted by a Levy distribution law at micro-scales. The relevance of a Levy distribution is 
made plausible by a simple model of nonlinear transfer. 

PACS numbers: 02.50.-r,02.50.Ey,64.60.aI 
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I. INTRODUCTION 



Many natural or man-made phenomena, such as turbu- 
lence flows, fluctuations in finance (stock market), seis- 
mic recording, internet traffic, climate change, etc., are 
characterized by randomness or stochasticity Ju [3, 0, 0, 

S i 0, i, S m [O, [13, [H, Q [H, [H, [M El- Anal- 
ysis of non-stationary stochastic processes referring to 
quantities which fluctuate widely and are uncertain has 
been a problem of fundamental interest for a long time. 
Over the past two decades, oil price has increased very 
sharply, rising from $25 per barrel in January 1986 to a 
peak of close to $ 122 per barrel in the last week of July 
2008. The effects of oil price fluctuations on the world 
economy are undeniable and particularly evident from 
the international reports. Oil price data as a time series 
exhibit complex patterns and seemingly appear to be a 
chaotic system. Indeed, the behavior of oil price fluc- 
tuations can be efficiently modelled by the Isingmodel 
which was proposed for stock-price fluctuations |19[ and 
by the the cascade model developed based on fractal con- 
cepts which was used for hydrodynamics and magneto- 
hydrodynamic turbulence [13, [SI, [l!] • Here, we employ 
the latter technique to characterize the statistical prop- 
erties of the oil-price time series, which sensitively dis- 
tinguish between self-similarity and multi-fractality in a 
time series. The model is based on two-point increments 
of oil price, yielding a comprehensive and scale-dependent 
characterization of the statistical properties of the system 
via an associated Probability Density Function(PDF). It 
is necessary to stress that the data series is represented 
by a finite number of records which do not constitute 
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a stationary process. The effect of non-stationarity on 
the detrended fluctuation analysis has been investigated 
in Ref. [1^. To eliminate the effect of sinusoidal trend, 
we apply the Fourier Detrended Fluctuation Analysis (F- 
DFA) . After the elimination of the trend we use the Mul- 
tifractal Detrended Fluctuation Analysis (MF-DFA) to 
analysis the data set. The MF-DFA methods are the 
modifled version of detrended fluctuation analysis (DFA) 
to detect multifractal properties of time series. The de- 
trended fluctuation analysis (DFA) method introduced 
by Peng et al. [23| has became a widely-used technique 
for the determination of (multi)fractal scaling properties 
and the detection of long-range correlations in noisy, non- 
stationary time series [2j, [2J] . It has successfully been 
applied to diverse fields such as DNA sequences [H, , 
heart rate dynamics [l^, [13 , neuron spiking , human 
gait [29! . long-time weather record s Iwi . Isd^ , cloud struc- 
ture jSlj ■ geology [3^ , ethnologv [33l|. economical time 
series 0, [SJ] , solid state physics |33 , sunspot time series 
[3^ . and cosmic microwave background radiation [T]].[l^. 

We propose a method for generating a stationary pro- 
cess analysis out of a non-stationary process. The fact is 
that stationary stochastic systems often show scaling in a 
statistical sense, coincident with non-Gaussian leptokur- 
tic (heavy- tailed) statistics. Importantly, identification 
of the associated scaling exponents implies the ability to 
interpret and estimate the behavior of the fiuctuations as 
well as the detection of long-range correlations. A self- 
similar Brownian walk with Gaussian PDFs, which has 
scaling exponent 1/2, is a good example of the process 
where shows uncorrelation on all temporal scales. We 
try to determine the scaling properties of the PDFs that 
are leptokurtic at micro-scales. The scaling exponents 
can be determined through the scaling behavior of the 
moments, usually characterized by computing structure 
functions. It is said that the fluctuations are self-similar 
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(monofractal) if scaling exponents of the moments exhibit 
a hnear power-law dependence. In contrast a nonlinear 
dependence infers to multifractal scaling, which is caused 
by intermittent small-scale structures of oil price fluctu- 
ations. A similar feature has been found in physical sys- 
tems for example, in velocity and magnetic fields of the 
solar wind [37] . l38l [39j as well as in magnetohydrodynamic 
turbulence studied via direct numerical simulations [2^ . 
Finally, distributed price changes are characterized by a 
stable Levy distribution in the central part of the distri- 
bution. 

The paper is structured as follows. In Section II we 
describe our data set. The MF-DFA method is briefly 
presented in Section III and shown that the scaling ex- 
ponent determined via the MF-DFA method are identical 
to those obtained by the standard multifractal formalism 
based on PDF analysis. In Section IV we employ a re- 
cently developed technique [13, [13, [1^ that sensitively 
distinguishes between self-similarity and multifractality 
in times series. By analyzing the temporal evolution of 
price dynamics, we demonstrate the strongly the non- 
Gaussian behavior of the returns of oil price and scale- 
dependent behavior of the PDFs. Also we explain the 
Hurst exponents analysis. The micro-scale PDFs resem- 
ble leptokurtic Levy distribution which will be discussed 
in Section V. Finally, in Section VI we will summarize all 
results discussed throughout this paper. 



II. THE DATA 
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FIG. 1: Semi-log plot of oil price series over the period 1986- 
2008. Inset: daily return of the oil price index. 



TABLE I: Mean, standard deviation, skewness, and kurtosis 
of the oil price returns. 



Mean 


Standard Deviation 


Skewness 


Kurtosis 


0.0124274 


0.729173 


-0.606749 


9.55225 



leptokurtic distribution), which confirms the existence of 
intermittency in the fluctuations. On the other hand 
the large value of kurtosis, k = 9.55225, with respect 
to Gaussian kurtosis {k — 3), show that the tails of the 
return distribution are fatter than the Gaussian ones. 



Over a twelve- year period, on average, the price of oil 
has increased from $25 per barrel in January 1986 to a 
peak of close to $ 122 per barrel in the last week of July 
2008 . Oil price as recorded in international markets 
[4Q] offer us an almost unique possibility to gain infor- 
mation on the stochastic dynamics state in a very large 
scale range, say 1 day up to 200 days. Our database 
consists of about 5695 daily price values which seem to 
provide a set of data points which will be sufficient to 
obtain the scaling properties the system. Fig[T] presents 
the daily fluctuations in oil price p{t) in the period 1986- 
2008. It is evident from the figure that the fluctuations 
do not constitute a stationary process, for instance one 
can show that the variance of the signal in some win- 
dow does not remain stable upon increasing the window 
size. Let us introduce the increments (or returns) Sp{t, t) 
defined by, Sp{t, t) = p{t + r) — p{t). The resulting se- 
ries for Sp{t) is shown in the inset graph of Fig[TJ It is 
straightforward to show, by measuring the average and 
the variance of Sp{t, r) in a moving window, that Sp{t, t) 
is stationary. Upon initiating the analysis of the distri- 
bution of oil price returns, the mean, standard deviation, 
skewness, and kurtosis of the return series are calculated 
(see Table |l|. It is easy to show that the skewness of 
a Gaussian function is zero so that the negative value 
of skewness, A = —0.606749, is a hallmark of departure 
of the PDFs from the Gaussian distribution (such as a 



III. THE MF-DFA ANALYSIS 

The MF-DFA method is a modified version of de- 
trended fluctuation analysis to detect multifractal prop- 
erties of a time series. Omitting unnecessary details, a 
brief summary of the method for calculating MF-DFA 
based on fractal concepts can be formulated in five steps. 
We take the price series pk with the size of N and follow 
the steps as follow: 

• step 1: Determine the " profile" 

i 

Y{t) = Y,[Pk'{p)], z = l,....,iV (1) 

fc=i 

where {p) is the mean of the series. Subtraction of the 
mean (jp) is not compulsory, since it would be eliminated 
by the detrending later in the third step. 

• step 2: Divide the profile Y{i) into Ns = int{N/s) 
nonoverlapping segments of equal lengths s. Since the 
length N of the series is often not a multiple of the con- 
sidered time scale s, a short part at the end of the profile 
may remain. In order not to disregard this part of the 
series , the same procedure is repeated starting from the 
opposite end. 

• step 3: Calculate the local trend for each of the 2Ns 
segments by a least-square fit of the series. Then deter- 
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FIG. 2: The fluctuation function F2{s) as a function of box 
size s for the returns of the oil index. 



mine the variance 

F^{s, i^)^- ^{r[(^ - 1). + ^] - (2) 

for each segment v, i> — 1, Ng and 
F'{s, v)^- Y.{Y[N -[y- Ns)s + ^] ^ (3) 



for V — Ng + 1, ...,2Ns. Here we use the Unear fitting 
polynomial in segment v. 

• step 4- Average over all segments to obtain the q-th 
order fluctuation function, defined as: 



2Ns 



9/2 



1/9 



(4) 



where, in general, the index variable q can take any real 
value except zero. We repeat steps 2, 3 and 4 on several 
timescales s. It is apparent that Fq{s) will increase with 
increasing s. 

• step 5: Determine the scaling behavior of the fluctu- 
ation functions by analyzing log-log plots of Fq{s) versus 
s. 

A power law relation between Fq{s) and s indicates 
the presence of scaling: Fq{s) ~ s"^'^\ In general, the 
exponent a(q) may depend on q. For stationary time se- 
ries such as fGn (fractional Gaussian noise), Y{i) in Eq. 
([T]), will be a fBm (fractional Brownian motion) signal, 
so, < a{q = 2) < 1.0 [4l|. As a result, the fluctua- 
tion function -^2(5) shows a scaling behavior, a(2), which 
is identical to the well-known Hurst exponent H. The 
Hurst exponent is called the scaling exponent or correla- 
tion exponent, and represents the correlation properties 
of the signal. If H — 0.5, there is no correlation and the 
signal is an uncorrelated signal, if iJ < 0.5, the signal is 



anticorrelated, if _ff > 0.5, there is positive correlation 
in the signal. We obtain the following estimate for the 
Hurst exponent, H — 0.52±0.02 as we can see in Fig. ©. 
Since H > 0.5 it is concluded that oil price returns show 
persistence, i.e, a certain correlation among consecutive 
increments. For s ^ 150 the empirical data deviate from 
the initial scaling behavior, as we can see in Fig. 
This indicates that oil price tends to loose its memory 
after a period of the order of 200 days or less. 



IV. STATISTICAL SELF-SIMILARITY 

A set of time series 6p{t, r) is obtained for each value of 
nonoverlapping time lag r. The return of the stochastic 
variable 6p{t, r) is said to be self-similar with parameter 
Q! (a >), if for any A 



dp{T) =^ X"6p{Xt). 



(5) 



The relation ^ is interpreted as an equality in law (EL), 
that is the two sides of the equation have the same sta- 
tistical properties. For the associated cumulative proba- 
bility distribution, it follows that 



piSpir) <p)^ piX'^SpiXr) < p), 



(6) 



for any real p. This implies for the probability density P 



P[Sp{T)] = X-"P,[X-''Sps], 



(7) 



introducing the master PDF Pg with 6ps — 5p{Xt). Ac- 
cording to Eq. (O, there is a family of PDFs that can 
be collapsed to a single curve P^, if a is independent 
of r. This is known as monoscaling in contrast to mul- 
tifractal scahng observed, e.g., for oil price returns at 
macro-scales. 

To characterize quantitatively the observed stochastic 
process, we measure the PDF P{5p) of the price fluctu- 
ations for T ranging from 1 to 200 days. The number of 
data in each set decreases from the maximum value of 
5695 (r = 1) to the minimum value of 5495 (r = 200). 
In Fig. ([3]) we show the four selected PDFs (normalized 
with the variance < 6p{t)'^ >^/^) for t — I, t — 20, 
T = 60, and r = 200 from the top ( left side) to the bot- 
tom (right side) respectively. The distributions lose their 
leptokurtic shape, as r increases. Due to the lack in cor- 
relation among distant fluctuations, the associated distri- 
butions become approximately Gaussian at macro-scale. 
The scaling behavior of the distribution at coarser time 
scales has two different regimes. At micro-scales (typi- 
cally shorter than 10 days), correlations between succes- 
sive price changes are dominant. This may be due to 
several reasons, such as oil pipe line damage, weather 
changes, or local variations in the internal (US) oil avail- 
ability. Interestingly, the PDFs at the micro-scales have 
the same leptokurtic shape, exhibit monoscaling, and do 
not change fundamentally, and resemble closely Levy dis- 
tributions, see Fig. On the other hand, at macro- 
scale (typically larger than one month) permanent crisis 
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FIG. 3: The PDFs of the normahzed oil price fluctuations for 
four different scales, (a): t — 1, (6): r — 20, (c): r = 60, and 
(d): T — 200. For comparison, the Gaussian distribution with 
the same variance is depicted (solid line) 
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FIG. 4: The PDFs of the oil price fluctuations on micro- 
scales. We can see that the shape of the PDFs do not change 
fundamentally as a sequent of the monofractality. 

in the Middle East and north African govern the price 
drift and corresponding to a Gaussian regime. This is 
coincide with a multifractal feature, as we can see in Fig. 

However the micro-and macro time scales regimes 
can be led to a linear and nonlinear scaling-dependence 
respectively. 

First, let us consider the scaling as defined by the struc- 
ture functions. The generalized structure function of or- 
der n are simply defined as 

/ + 00 
\5prP{Sp,T)diSp). (8) 
-oo 

The analysis which follows is also valid for the moments; 
however, structure functions are typically calculated for 
a data series. The arguments do not apply to structure 



FIG. 5: The PDFs of the oil price fluctuations on macro 
scales. The PDFs become Gaussian as the time scales is in- 
creased. This is coincide with multifractality feature. 

functions of odd order, which not only may have negative 
coefficients, but could in fact even change the sign of the 
scaling range. The proof will, however, remain valid for 
odd orders if the structure functions are defined with the 
absolute value of the returns. Using the relation ((7]) we 
obtain 

5"(t;±oo) = X<-S:{Sp,;±oo), (9) 

where the linear function C„ — an refer to the statis- 
tical self-similarity, monoscaling case. On the contrary, 
in some cases, one may observe multifractality scaling, 
in the sense that a nonlinear dependence is observed 
on n where Cn = na{n) is a convex function of n and 
Cn+i > Cnyn. This deviation from strict self-similarity 
over all time scale t, also termed multifractal scaling, is 
caused by the intermittent micro-scale structure of tur- 
bulence. 

To test if the above-mentioned interesting observations 
in oil price are a phenomenon related to inherent proper- 
ties of stochastic processes, structure functions S"{t) for 
different r, given by Eq. ([5]), are computed. A difficulty 
that can arise in the experimental determination of the 
C„ is that for a finite length times series, the integral 
Eq. ^ is not sampled over the range (— oo; -f oo), rather 
the outlying measured values of y determine the limit, 
[—y;+y]. Fig. ([6]) shows some selected 5'"(t), according 
to Eq. ([5]). The slope of the curve gives scaling exponent 
C„ which can be obtained by fitting a straight line to 
a log-log plot in the interval / e [1,100]. Because of 
increasing of the statistical errors at the higher orders, 
such a fitting becomes rather arbitrary. In Fig. ([7]) 
we report the scaling exponents extracted from the 
structure functions. The behavior of ^„ against n shows 
that scaling exponents have nonlinear behavior at all 
scales , say, they are different from the usual linear an 
law. 
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FIG. 6: The structure functions are depicted, as computed 
from Eq. (|8}. In order to obtain the scaling exponents, 
we take the logarithmic slope of linear least-square fits (solid 
line) . 
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FIG. 7: The scaling exponents of the structure functions 
are depicted versus the corresponding order n. The nonlinear 
behavior is obvious in the figure. 



To apply the rescaling procedure given by Eq. ([7]) the 
exponent a is extracted from the underlying data by a 
independent technique, (2^. Ibt}. The standard deviation 
which is defined by the root of the second-order structure 
function, CT(r) = [S'^(t)]^/^ and has the minimum of sta- 
tistical error, exhibits power-law behavior with respect 
to the increment distance, a^r) ~ r" as depicted in Fig. 

A linear least-square fit is carried out to obtain a. 
The characteristic exponent deduced in this way is a = 
0.52 ±0.031 which is in good agreement with the value of 
iJ obtained in Section. jlTD). Fig. ([9]) and Fig. (fTO)) show 
the rescaled sets of PDFs on the micro- and macro-scales 
respectively. Evidently the PDFs at micro-scales are self- 
similar and collapse with weak scattering on the master 



FIG. 8; Standard deviation of oil price returns within the 
desired range. 



PDF, Ps, when using the characteristic exponents given 
above. The corresponding increment distances r arc all 
shorter than 10 days. We may model this PDFs by a 
Levy distribution, which thus turns out to be a success- 
ful fit to the distribution of oil price fluctuations. On 
the other hand, at macro-scale the PDFs do not show 
a self-similar behavior and rather constitute a multicas- 
cade process. This occurs at scales larger than the one 
month. Because of the resulting multifractal scaling of 
the PDFs it is evident that they can not collapse onto a 
single curve, see Fig. ([T0|) . 



V. LEVY DISTRIBUTION MODEL 

Levy-stable laws are a rich class of probability distri- 
butions that comprise fat tails and have many intrigu- 
ing mathematical properties. They have been proposed 
as models for many types of physical and economic sys- 
tems. There are several reasons for using Levy-stable 
laws to describe complex systems. First of all, in some 
cases there are solid theoretical reasons for expecting a 
non-Gaussian Levy stable model, which can be a good 
fitting to experimental and numerical results. The sec- 
ond reason is the Generalized Central Limit Theorem 
which states that the only possible non-trivial limit of 
normalized sums of independent identically distributed 
terms is Levy-stable. (Recall that the classical Central 
Limit Theorem states that the limit of normalized sums 
of independent identically distributed terms with finite 
variance is Gaussian.) The third argument for modeling 
with Levy-stable distributions is empirical; many large 
data sets exhibit fat tails (or heavy tails); for a review 
see [12, . Such data sets were described by a Casting 
model based on the lognormal ansatz (in terms of 
the variance of the Gaussian distribution). To con- 
firm that model it is convenient to summarize the basic 
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and the symmetric Levy distribution becomes 
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FIG. 9: Rescaled PDFs of oil price fluctuations in micro-scales 
. The Levy law is represented by the dashed line. 
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FIG. 10: Rescaled PDFs of oil price fluctuations in a macro- 
scale. This shows that the collapse of the PDF's to a single 
curve is broken over the scales. 



features of the Levy stable distribution. 

A Levy process is a time-dependent or position- 
dependent process that at an infinitesimal interval has 
the Levy distribution of the process variable. The char- 
acteristic function of the Levy process is 



K^{q,s) = exp(-cs | q |^) , 



(10) 



where s can be a characteristic time or space scale. If /i = 
2 the Levy collapses to the Gaussian distribution. If = 
1 the Levy becomes a Cauchy distribution. The original 
Levy process is given by its inverse Fourier transform, i.e. 



Pf_,{x,s) = J dqe 



iqx—cs\q\^ 



(11) 



1 



exp{—-fAsq'^)cos{qSxAs)dq, (12) 



where the increment is Sx = Xs - 
2, and 7 > is a scale factor, 
probability leads to 



Xs-As] here, < /i < 
The maximum event 



r(i//i) 

7r^(7As)i/A' 



(13) 



The exponent fj, of the best fits is constant at the micro- 
scales range and amounts approximately to /i ^ 1.92 
which is fi = 1/a. Natural phenomena also investigated 
where similar findings have been reported include, for 
example, financial systems e.g. the Tehran price stock 
market, where /i ~ 1.36 [1], and also physical systems 
such as the solar wind, where fi ~ 3.3 [44{ . From Fig. 
^ we conclude that a central section of Levy distribu- 
tions describe very well the dynamics of the PDFs of oil 
price fluctuations at micro-scales. One can see that the 
rescaled PDFs are definitely non-Gaussian. 



VI. SUMMARY 



In this paper, we have presented a statistical analysis 
of oil price in the United States for the period of 1986 
to 2008. We have apphed a generic MF-DFA method to 
extract scaling exponent of the fluctuation functions, in 
particular, relying on the second order ^2(5) which was 
used in the rescaling procedure. However, the simple 
scaling properties that we have found via an analysis of 
the PDFs, allow us to detect mono(multi) fractality fea- 
ture over all time scales. The presence of intermittency 
in oil price fluctuations is manifested by the leptokurtic 
nature of the PDFs which show increased probability of 
large fluctuations compared to that of the Gaussian dis- 
tribution. Fluctuations on the macro temporal scales, 
T > 10 day, converge toward a Gaussian distribution and 
are an uncorrelated signal. The reason may be due to 
unstable conditions in the Middle East or OPEC's de- 
cisions oil production. The critical macro scale which 
was obtained is different for some financial and physical 
systems; see for example in Refs. [1, [l^l- We have also 
obtained a good collapse onto a single curve for r < 10, 
according to the rescaling procedure O. The proximity 
of the PDFs to Levy distributions is made plausible by a 
simple model mimicking nonlinear spectral transfer. 
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